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In  provlous  theoretical  treatments  of  the  stability  of  the  compressible 

i  . 

•••  laminar  boundary  layer  the  effect  of  the  temperature  fluctuations  on  the 
J  "viscous"  (rapidly-varying)  disturbances  is  either  Ignored  (Loes-Lln),  or 
S  Is  accounted  for  Incompletely  (Dunn-Lln).  A  thorough  reexamination  of 
this  problem  shows  thar  temperature  fluctuations  havo  a  profound  influence 
on  both  the  "invtscld"  (slowly-varying)  and  viscous  disturbances  above  a 
Mach  number  of  about  2,  0,  -The  present  analysis  includes  the  effect  of 
j  temperature  fluctuations  on  tho  viscosity  and  thormal  conductivity,  and  also 
Introduces  the  viscous  dissipation  term  that  was  dropped .ln-tha~eaxllar 
theoretical-  treatments. 

Some  important  results  of  tho  present  study  are:  (1)  Instead  of 
j  being  nearly  constant  across  the  boundary  layer  the  amplitude  of  the 
j  invlscid  pressure  fluctuations  decreases  markedly  with  distance  outward 
j  from  the  plato  surface  for  Mach  numbers  greater  than  3.  This  behavior 
I  means  that  the  Reynolds  shear  stress  near  the  critical  layer  Is  greatly 
reducedj'v(2)  at  Mach  numbers  less  than  about- 2  dissipation  effects  are 
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minor,  but  they  beconio  oxtromcly  Important  at  higher  Mach  numbers  and 
tend  to  becomo  the  dominant  effect  for  hypersonic  Mach  number  if* 

(if  the  rate  of  conversion  of  enorgy  from  tho  mean  flow  to  the  disturbance 
flow  through  the  action  of  viscosity  at  tho  wall  Increases  with  Mach  number^ 
(4)  the  minimum  critical  Reynolds  number  for  Insulated  flat  plate  boundary 
layers  deci  cases  In  tho  range  0  *  M  »  3,  and  then  rises  vary  sharply  for 
hypersonic  Mach  numbei 

Numerical  exampNrrillustrattng  tho  effects  of  compressibility 
(including  neutral  stability  characteristics)  are  obtained  at  Mach  numbers 
of  2.  2  and  5. 6.  Tho  calculated  neutral  stability  diagrams  are  compared 
with  the  experimental  results  of  Laufor  and  Vrebalovlch  at  M  ■  2.  2,  and 
of  Dcmetrladoe  at  M  ■  5.  3. 
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quantity  of  the  total  flow 

mean  or  steady  component  of  flow  quantity 

fluctuating  component  of  flow  quantity 
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gas  constant 
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ratio  of  specific  heats 
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small  parameter  (aRe)  * 
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boundary  layer  momentum  thickness 
wave  length  of  disturbance 
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A  bar  over  a  quantity  Indicate*  mean  value. 

Prime*  generally  denote  differentiation  with  respect  to  y.  The  few  instances 
where  prime*  denote  a  fluctuating  quantity  should  not  cause  any  confusion. 

1.  Introduction 

Experimental  result*  obtained  by  Laufer**  and  Vrebalovlch  at 
■  2.  2  furnish  definite  proof  of  the  existence  of  Tollmlen  waves  In  the 
supersonic  laminar  boundary  layer  on  an  Insulated  flat  plate.  The  general 
shape  of  the  stability  diagram  tn  the  ag  -  Reg  plane  does  not  differ  much 
from  that  found  at  low  speeds;  however,  the  minimum  critical  Reynolds 
number  and  the  wave  numbers  and  amplification  rates  of  the  self-excited 
disturbances  arc  much  lower.  Tollmlen  waves  were  also  found  by  Demetrlades* 
in  the  laminar  boundary  layer  of  an  Insulated  flat  plate  at  M  3  5.  8.  Here  the 
minimum  critical  Reynolds  number  seems  to  be  about  one  order  of  magnitude 
larger,  while  the  characteristic  wave  numbers  and  amplification  rates  are 
even  lower  than  atM  :  2.2, 

These  experimental  studies  stimulated  a  reexamination  of  the  whole 

theoretical  basis  of  laminar  stability  at  high  Mach  numbers.  In  the  earlier 
3  4 

treatments  of  Lin  *  and  one  of  the  present  authors,  it  is  tacitly  assumed 
that  the  critical  layer  (wave  speed  =  flow  speed)  lies  close  to  the  surface, 
so  that  the  flow  velocity  with  respect  to  the  wave  is  small  in  the  region 
between  the  surface  and  this  layer.  In  that  case  the  rate  of  change  of 
density  (or  temperature)  following  a  fluid  particle  is  small  for  the  "viscous" 


+  Superscripts  refer  to  references  listed  at  the  end  of  the  paper. 


(rapidly-varying)  disturbances,  and  thc»o  disturbances  are  treated  as 
ossontlally  Incompressible.  Dunn5  and  Lin  pointed  out  that  this  conclusion 
Is  valid  for  subsonic  or  slightly  supersonic  speeds,  but  not  for  higher 
supersonic  speeds,  because  the  critical  layer  is  no  longer  "close"  to  the 
surface.  They  obtained  the  "first-order"  offoct  of  temperature  fluctuations 
on  the  viscous  disturbances  through  the  continuity  equation,  and  they 
discussed  the  Importance  of  the  thermal  boundary  condition.  However, 
only  the  so-called  "leading  terms"  are  retained  In  the  energy  and  momentum 
equations.  No  numerical  calculations  utilising  this  method  are  available. 

At  high  Mach  numbers  the  temperature  fluctuations  are  dominant, 

and  viscous  dissipation  must  become  important.  In  the  present  study  we 

want  to  retain  these  effects,  so  far  as  possible.  Even  for  the  "invlscld" 
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(slowly-varying)  disturbances  previous  treatments  have  employed 
approximate  methods  valid  only  up  to  moderate  supersonic  speeds,  and  for 
wave  velocities  not  too  close  to  free- stream  speed.  New  methods  must 
now  bo  developed  to  cope  with  high  supersonic  and  hypersonic  speeds. 

Only  the  simplest  gas  Is  considered  here,  namely  one  with  constant 
specific  heats,  constant  Prandtl  number  and  constant  chemical  composition. 
Of  course  the  temperature  dependence  of  viscosity  and  thermal  conductivity 
is  taken  Into  account.  The  laminar  boundary  layer  la  idealized  as  a  plane, 
parallel  flow,  and  each  total  flow  quantity  Is  regarded  as  being  composed  of 
a  mean  component  which  dopands  only  on  the  distanco  normal  to  the  plate 
surface,  and  a  fluctuating  component  of  infinitesimal  magnitude.  Thus 
Q(*.y.  t)  ■  CF  (y)  ♦  Q'  (x,  y,  t) 

With  these  assumptions  the  coefficients  of  the  linear  partial  differential 
disturbance  equations  are  independent  of  both  x  and  t,  so  that  a  disturbance 
of  the  form 

Q'  (x,  y.t)  a  q(y)  eio^X  '  Ct^ 

is  suggested.  Here  q(y)  is  a  complex  amplitude  and  c  is  the  comptex  phase 
velocity.  The  real  part  of  c  =  c  is  the  dimensionless  velocity  of  wave 
propagation  parallel  to  the  surface,  while  ac^  gives  the  amplification  (or 
damping)  rate. 

The  total  flow  must  satisfy  the  same  boundary  conditions  on  the 
velocity  as  the  original  steady  mean  flow,  so  that  the  longitudinal  and  normal 
velocity  fluctuations  vanish  at  the  surface,  i.  e.  ,  fw  a  0  and  a  0,  In 
general  the  thermal  boundary  condition  at  the  surface  states  that  the 
instantaneous  temperature  and  heat  transfer  rate  are  continuous  across  the 
solid-gas  interface.  However,  most  surface  materials  are  so  highly 


conductive  compared  to  gases  that  the  temporaturo  fluctuation*  at  the  surface 
are  almost  completely  suppressed  at  tho  frequencies  of  interest  for  laminar 
stability.  Therefore  we  take  0^  ■  0.  In  Section  2  wo  show  thAt  for  disturbances 
propagating  at  subsonic  velocities  relative  to  tho  free  stream  alt  disturbance 
amplitudes  vanish  far  from  the  plate  surface,  l.  e. ,  q(y)  — *0  as  y— 

The  purpose  of  the  present  paper  Is  to  bring  out  the  physical 
mechanisms  and  main  theoretical  problems  of  laminar  boundary  layer 
stability  at  high  Mach  numbers.  Details  of  the  mathematical  treatment 
(Including  the  general  thermal  boundary  condition  a  O  '  I  b  •  0)  are 
contained  in  Reference  6.  In  Section  2  we  detineate  the  principal  disturbance 
flow  parameters  and  regions  In  coordinates  fixed  in  the  wave.  Section  3 
is  concerned  with  the  mechanisms  of  production  of  disturbance  energy  at 
high  speeds,  and  the  conclusions  to  be  drawn  from  the  energy  balance  for 
neutral  disturbances  (c^  ■  0).  In  Section  4  some  aspects  of  the  eigenvalue 
problem  are  discussed,  with  particular  emphasis  on  the  role  of  the 
temperature  fluctuations,  and  some  numerical  examples  are  presented 
for  an  Insulated  flat  plate  In  air  at  »  2.  2  and  5.  6.  Finally,  in  Section  5 
we  summarise  our  conclusions  and  the  present  state  of  our  knowledge 
(and  ignorance)  of  laminar  stability  at  high  Mach  numbers. 

2.  Flow  Regions  and  Parameters  in  Coordinates  Fixed  in  the  Wave 

To  an  observer  riding  with  the  wave  the  entire  flow  field  Is  steady. 

The  uniform  free  stream  is  moving  to  the  right  with  the  velocity  1  -  c  ,  while 
the  plate  surface  la  moving  to  the  left  with  the  velocity  c  .  (See  sketch  on 
page  7.)  The  present  discussion  is  concerned  only  with  disturbances 
propagating  at  subsonic  velocity  relative  to  the  free  stream.  In  other  words, 
the  relative  Mach  number  (MREL)fl  ■  M#  (1  -  cf)  a  1,  or  cr  »  l  -  (l/M^). 

In  such  a  steady,  subsonic  flow  of  unlimited  extent  it  is  well-known  that  all 
small  disturbances  die  out  with  distance  liko  e”^  as  y  — foe  ,  where 

p  *  e  /  1  -  (MR£L)e2  *  a  fl  -  Me2(l-cr)2  . 

The  restriction  to  subsonic  relative  motion  evidently  does  not  apply 
to  the  plate  surface.  In  wave  coordinates  the  plate  Mach  number  is  given  by 
(M^ee)w  »  Mfl  c^  /^Tw<  For  an  insulated  surface  and  Pr  »  1,  one  finds 

(Mfl  c  )/(  ^  Tw  )  ^  1  when  ^  2.  2  (approximately,  for  Y  a  1.4).  Thus 


for  Mc  <  2.  2  the  flow  ia  everywhere  subsonic  with  respect  to  the  wave,  but 
for  Mc  >  2.  2  a  supersonic  flow  region  exists  near  the  plate  surface.  In  fact, 
when  >  >  I,  c,  — 1  ,  and  (Me  cr)/(  fr^  )  .  In  this  limiting  case 

the  sonic  line  defined  by  the  relation  Me2  (w  -  c)2  a  T  occurs  at  w  a  2/3*. 

The  existence  of  an  extensive  supersonic  flow  region  for  sufficiently 
high  Me  has  a  profound  effect  on  the  amplitude  distribution  for  the  "Inviscld" 
pressure  fluctuation  between  the  plate  surface  and  the  critical  layer  |w  >  c). 
At  the  plate  surface  =  0*  therefore.  In  wave  coordinates  the 

amplitude  of  the  pressure  fluctuations  near  the  plate  surface  satisfies  the 
Prandtl-Glauert  equation  for  a  wavy  disturbance  in  a  steady,  uniform  flow, 
namely, 

*M  *  °  [MREL  "  *]  *  a  0  •  (1) 

From  Eq.  (1)  one  sees  that  the  amplitude  of  the  pressure  fluctuation 
decreases  Initially  with  distance  away  from  the  plate  surface  if  M  >1 


*  This  sonic  line  does  not  introduce  any  real  singularities  into  the 
linearized  disturbance  equations  for  a  shear  flow'*  . 


and  increase*  if  M,JPI  <  1.  If  the  y-coordlnnte  la  normalised  by  mean*  of 
RE.v  Y 

the  liowarth-DorodnlUyn  transformation  Y  ■  \ 


then  Eq.  (1)  become** 


d2* 


♦  a2  T 


w  frMREL  ^w  ‘  *1 


(2) 


In  other  word*  the  proper  parameter  I*  not  a2»  but  n2  T. 


.  2  ■  ■-  Z  i  *■"  • 

w  M8  for 

an  insulated  plate.  At  subsonic  or  slightly  supersonic  free  stream  speed* 
a2  T  2  *  6(0*1  •  and  the  increase  in  pressure  amplitude  between  the  plate 
surface  and  the  critical  layer  is  small.  But  at  high  supersonic  and  hypersonic 
speeds  «2T  2  is  no  longer  small,  and  ®o  that  the  decrease  in 

pressure  fluctuatlpn  amplitude  outward  to  the  critical  layer  is  substantial. 

This  phenomenon  Is  not  properly  accounted  for  In  previous  theoretical 
treatments'"^,  because  the  tnvlscid  solutions  are  obtained  by  employing 
series  expansions  in  powers  of  a2.  Since  a2  is  supposed  to  be  small,  it 
is  tacitly  assumed  that  |»c/»w  |2  ■  Of  1 ).  Tho  fact  that  |*c/ww  |2  <  <  1  at 
high  Me  has  a  strong  influence  on  the  Reynolds  stress  increment  (or  decrement) 
at  the  critical  layor,  and  therefore  on  the  energy  balanco  for  a  neutral 
disturbance  (Section  3). 

The  distinction  batween  subsonic  and  supersonic  flow  regions  in 
coordinates  fixed  In  the  wave  makes  sense  only  if  the  disturbances  aro 
largely  "tnvtsc Id'*  (slowly-varying),  as  proposed  originally  by  Prandtl. 

Viscosity  and  conductivity  aro  important  In  two  regions:  ( 1 )  at  the  plate 
surface,  where  viscous  solutions  must  be  added  to  the  tnvlscid  solutions  in 
order  to  satisfy  the  boundary  conditions  fw  ■  0,  ®w  *  0»  (2)  at  the 
"critical  layer"  (w  ■  c),  where  the  longitudinal  transport  of  vortlclty  and 
heat  energy  with  respect  to  the  wavo  vanishes,  and  the  vertical  transport  of 
these  quantities  must  be  balanced  by  viscous  diffusion  and  heat  conduction. 
Clearly  the  proper  local  Reynolds  number  in  coordinates  fixed  in  tho  wave  is 


aRe 


aRe  c 


|  (w  -  c  )|  .  At  tho  plate  surface  the  parameter  (  — - -  )  *  Is  a 

r  vw 


measure  of  the  diffusion  distance  for  vortlclty  during  one  period,  *+  so 
this  inner  boundary  layer  is  thin  when  oRe  >  >  i.  Similarly,  the  parameter 


*  Since  (0INV)W  =  0. 


*w=°- 


*+  The  parameter  {  a^<~  ^  ,  )"^  measures  the  corresponding  diffusion 

*  w 

distance  for  heat  energy,  where  <r  is  the  Prandtl  number* 


\w 


INV'w 


■  1*  *1  1‘inv' 


1.  e.,  tho  "invlscld"  temperature  fluctuation*  normalized  by  th«  fro*  atream 
temporatura  *ro  of  order  M  Because  of  tho  boundary  condition 

H 

*°lNV*w  *  ^Vw  "  0  the  "viscous"  temperature  fluctuation*  are  alao  of 

2  .  ■  '  ■  .  '  m  . : _ ■  1  _  # _ > -  .  . .  i 


Now  the  gradient  of  mean  temperature  1*  of  order  M#  and 

.  i  .  .  .  .  i  .  fti„i_..  _ .  .  . _  » .  1 


order  M4 

*o  1*  the  whole  level  of  mean  temperature  in  the  boundary  layer  when  M#*  >  >  1 
In  other  word*  the  free  stream  static  temperature  Is  no  longer  relevant,  and 
must  be  replaced  by  some  representative  temperature  T#  » 


For  example,  the  new  parameter  that  orders  the  viscous  terms  In  the 

1  Re 

equations  of  motion  Is  given  by  —  ,  where  Re___,  ■  ■  ■»  ■■  I 

<aReREF)*  REF  yREF 

therefore,  £  «  £  ,  where  £  Is  tho  old  parameter  (aRo)"^.  Fora 

linear  viscosity  temperature  relation  £-«M#2  .  Terms  in  the  equations  of 
motion  arising  from  tho  normal  gradients  of  mean  or  fluctuating  viscosity 
and  conductivity,  or  from  the  viscous  dissipation,  were  regarded  as  of  order 
£  In  References  3,  4,  and  5  compared  with  tho  "leading"  viscous  terms, 
but  they  are  actually  of  order  £  .  At  high  supersonic  or  hypersonic  free 
stream  speeds  these  terms  aro  likely  to  be  of  tho  same  order  as  the  so-called 
"leading  terms"  (Section  4). 


3,  Energy  Balance  and  Reynolds  Stross  Distribution  for  a  Neutral  Disturbance 

For  a  neutral  disturbance  the  net  rate  of  transfer  of  energy  per  cycle 
to  the  disturbance  by  the  action  of  tho  Reynolds  stross  Z  =  -  p  u*  v'  must  be 
balanced  exactly  by  tho  rate  of  dissipation  of  disturbance  energy  per  cycle. 

As  expected,  the  dissipation  term  is  linear  in  the  viscosity  and  is  therefore 
of  order  l/aRe.  It  is  somewhat  less  obvious  that  the  Reynolds  stress  also 
depends  on  viscosity  because  of  tho  phase  shift  in  disturbance  velocities 
Introduced  near  the  plate  surface.  Tho  Reynolds  stress  distribution  across 
the  boundary  layer  for  aRo  >  >  1  (Section  2)  is  sketched  on  page  11.  [  The 

rate  of  transfer  of  energy  to  the  disturbance  is  given  by  j  £?{ du/dy)  dy  dx  ] 

The  Reynolds  stress  rises  rapidly  with  distance  away  from  the  plate  surface 

[777- 

in  a  zone  of  thickness  $w  "i  ,  and  then  remains  practically 

constant  in  the  nearly-inviscid  region  between  the  surface  and  the  critical 


layer.  Since  t  — ♦»  0  far  from  the  surface,  V  must  bo  cancelled  by  an  equal 
and  opposite  jump  in  Reynolds  stress  at  th »  critical  layer.  When 

~~  ( l  •  c)  >  >  1  t  this  Jump  can  be  calculated  from  the  Inviscid  solutions 

(Soction  2);  viscosity  and  conductivity  merely  smooth  out  tho  transition.  Thus 

fw  r«*‘° 


»  yc  -  o 

Let  us  now  examine  each  of  theso  regions  of  Reynolds  stress  production 
(or  destruction)  separately. 


(a)  At  the  plate  surface 


For  a  fictitious  inviscid,  non-conducting  gas, only  the  inviscid  solutions 
apply.  All  disturbances  vanish  far  from  the  surface,  and  the  normal  velocity 
fluctuation  vanishes  also  at  the  plate  surface.  However,  the  inviscid  tem¬ 
perature  and  longitudinal  velocity  fluctuations  generally  take  on  non-zero 

values  at  the  surface.  Now,  for  a  real  gas  f  a  0  and  0  =  0,  no  matter 

w  w 

how  small  the  viscosity  and  thermal  conductivity.  Thus  one  must  add  viscous 
solutions  determined  by  the  conditions  (fy)w  a  -  (fjNV)w  and  (9y)w  »  -  (°INy)wS 
®v  — 0  as  y  — >  co.  (See  sketch  on  page  12. )  By  considering  only  the 
leading  viscous  terms  in  the  momentum  and  energy  equations  in  a  thin 
layer  near  the  plate  surface,  and  recognizing  that  0  and  w^  0  there,  one 
finds  the  approximate  viscous  solutions 


0-U(y  7S76w) 


where  6, 


But  these  additional  temperature  and 

longitudinal  velocity  fluctuations  induce  an  increment  in  normal  volocity 
fluctuation  across  the  boundary  layer.  According  to  the  equation  of 
continuity  and  the  equation  of  state 

8v*  f  3u*  1  dp*  1  dT*1 


-  -  i(v  -  (ic/TJ  0V  . 

By  utilizing  Eqs.  (S)  and  (9)  in  Eq.  (11),  employing  Eq.  (6) 
recognizing  that  3  0 


one  obtains 


J  <V  dy  -  ttfv)w  -  ■  (0V)W 

i<i  ♦  o  6W 


(12) 


^INV^w  ^  *  Kflr  )  » 


whore 


U-  1)M„2C2  , 

K  *  - “  J^MREL^w 


(»3) 


But  ■  y 

$w  *  <<Uw  +  ^INV^w  *  0 

■o  that 

^INV*w  *  *  «*v>w  » 

whore  (0y)w  1*  given  by  Eq.  (12),  and  the  tnvlscld  solutions  must  now  be 
altered  slightly  to  satisfy  the  boundary  condition  at  the  wall.  The  situation 
hero  is  somewhat  analogous  to  the  strcamllno  displacement  effort  produced 
in  tho  external  invlacld  flow  by  the  mean  boundary  layer. 

It  is  prcciacly  tho  fact  that  (0jjjy)w  /  0  for  oRo  l*rgo  but  not  Inflnito 
that  gives  rise  to  tho  Reynolds  stress  associated  with  the  lnviscld  solutions. 
Now  X  *  *  p  u'  v*  •  so  that 

V  -  -  (p/2)  U  Ri  (f  0  ♦)  .  (15) 

By  employing  Eqs.  (12)  -  (15),  one  finds 

r,NV  |Sw*  ’  I  'W'w  I  0  ♦«*-!)  .  (16) 

f Note  that  is  >  0  and  of  order  - 5 —  .  As  the  Mach  number  M 

L  ■  -/okz'/ :r "  - ^ 

increases  K  increases  ^Eq.  (13)]  ,  and  the  Reynolds  stress  near  the  plate 

surface  increases. 


(b)  Across  the  critical  layer  (w  =  c) 


The  magnitude  of  the  jump  in  Reynolds  stress  across  the  critical 
layer  is  given  by  the  behavior  of  the  inviscid  solutions  there^: 


AT 


=  -  (a/2)(c/Tw  Ww')  (fTMV) 


INV  w 


V  (c) 
o'  ' 


where 


y  -  0 
3z 

v  (c)  = 

o'  ' 

w  f 

>  TT  W 

T“ 

■>  w  ■ 

V  12 

Tf 

W 


(it) 


w  (wc  )' 


CTc  [ (d/dy)(w>/T)  ] 


w  =  c 
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<5)  connoctiiig  |f1Ny)^  nnil  ww 


limit  a  Re 


riNV  <V 


already  boon  utilized.)  In  the 
0  [Eq.  (16)]  ■,  and  the  Reynolds  stress 


vanishes  everywhere.  Thus  the  necessary  (nnd  sufficient)  condition  for  the 
existence  of  a  neutral  invlscld  disturbance  Is  that  (d/dy)(w'/x)  ■  0  for  some 


>  , 
w  ■  1  - 


)  (Reference  3).  The  value  of  w  ■  cg  for  Which  this  Condition 
is  satisfied  is  shown  in  Figure  1.  When  c  >  c  *  (d/dy)(w*' 


<  0  and 


yc  +  0 

«l 

#c-° 


<  0  ,  so  that  a  neutral  disturbance  can  in  fact  exist  when  aRa 


is  largo  but  finite.  Evidently  the  magnitude  of  this  jump  in  T  depends 
critically  on  the  ratio  |  «c/»w  P  • 

In  Section  2  it  was  shown  that  tho  amplitude  of  tho  pressure 
fluctuation  decreases  with  distance  away  from  the  plate  surface  at  high 
supersonic  spoeds.  A  quantitative  estimate  of  this  dccroaeo  (a  obtained 
from  tho  solutions  of  the  invlscid  equations  for  a  neutral  disturbance 
By  employing  Eqs.  (3)  and  (10)  ono  obtains 

0 


(c  -  c8). 


-  i  f 


UyjEL*  (i8) 

and  the  momentum  equations  parallel  and  normal  to  the  plate  surface  ylold 
ip(w-c)f+pw>0  -  -  —t  (19) 


r< 

1  a  2  p  (w  -  c)  0  m  -  (w'/y*M  2) 


(20) 


By  eliminating  f  and  0  ono  obtains  a  single  differential  equation  for  ir.  It 
turns  out  to  bo  more  convenient  to  introduce  the  Riccatl  transformation 

G  »  (w'/a2  w)  a  (l/a2)(d/dy)  log  n  , 
and  the  differential  equation  for  G  is  (Reference  6) 

M 

U 

- T 


l  - 


2  2  1 
,  (w  -  cf 


+  ( 


2w' 

w-c 


T*  2  2 

-ijr—  )  G  -  a  G 


(21) 


At  the  plate  surface  G  =  0  when  aRe  — >- «  and  we  obtain  Eq.  (1)  of  Section  2 
In  Figure  2  the  results  of  a  number  of  numerical  integrations  of 
Eq.  (21)  are  shown  for  an  insulated  plate  surfaco,  baaed  on  the  mean  flow 
profiles  calculated  by  Mack^.  The  location  of  the  critical  layer  coincides 
with  the  most  outward  zero  of  G.  The  behavior  of  the  slope  at  y  =  0  follows 
the  prediction  made  in  Section  2.  Since  tho  area  under  each  curve  out  to 
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any  distance  from  the  surface  In  proportional  to  tho  logarithm  of  tho 
presnure  fluctuation  nniplituilo  at  that  location,  it  is  clear  that  (»  /ir  )2 
decreases  markedly  with  increasing  Mach  numbor  at  high  MQ.  In  Plguro  3 
tho  ratio  »c/ww  in  plotted  an  a  function  of  MQ.  Below  M0  ^ 2.  5,  (irc/ww)  ■0(1). 
The  integration*  also  yield  tho  eigenvalues  (a.)  ,  and  those  values  are  shown 

In  Figure  4.  The  valuo  of  (og)M  is  very  small  when  MQ  <  <  1,  incrcasos  in 
the  subsonic  range,  and  (after  a  dip  near  Mfl  ■  1. 7),  reaches  a  maximum  at 
about  M#  »  5.  0.  Boyond  this  Mach  number  («q)#  approaches  the  asymptotic 
behavior  (n g)#  l/M#2  .  • 

By  referring  to  Eq.  (17)  and  Figure  3,  one  sees  that  the  stabilising 
(or  destabilizing)  action  of  tho  jump  in  Reynolds  stress  at  tho  critical  layer 
decreases  rapidly  at  high  Mach  number.  According  to  Eos.  (7),  (16),  and 
(17)  the  condition  for  a  neutral  disturbance  is 


This  relation  is  strictly  applicable  only  when  aRe  >  >  1,  I.  e. ,  along  the 
•'upper  branch"  of  the  curvo  of  neutral  stabMtty  in  the  o  -  Re  piano.  When 
M0  <  I  "  »  anti  Ihc  principal  effects  of  increasing  Mach 

numbor  are  to  shift  cf  to  higher  values  (cr  ■  1  -  (1/Mj  ,  and  to  Incroaso 
the  rate  of  production  of  disturbance  energy  near  tho  surfaco.  As  indicated 
bY  factor  (1  ♦  K  <J“  )  in  Eq.  (22),  this  last  effect  shifts  the  eigenvalues 
along  tho  upper  branch  to  still  larger  values  of  aRe,  and  is  essentially 
destabilizing.  ■  ■■■ '  -  : , 

:  At  hypersonic  speeds,  on  the  other  hand,  the  behavior  of  the 

quantity  vq(c)  J  *c/ww  j  is  the  dominant  factor.  The  variation  of  this 
quantity  with  cy  is  indicated  schematically  in  the  accompanying  sketch. 

Of  course  vQ(cg)  a  0,  and  the  product  v  (c)  j  sc/jtw  |2  at  first  increases 
with  increasing  cg  >  cg  .  But | vQ/ irw  |  <  <  l  (Figure  3),  and,  moreover, 

I  ffc/,rw  |2  dccrea sc 3  with  increasing  cr.  Thus  this  product  is  always  very 
small  numerically,  reaches  a  maximum  at  some  value  of  c  >  c  ,  and 
then  decreases  in  value  again  as  cr  — >  1.  According  to  Eq.  (22)  this 
behavior  means  that  aRe  decreases  from  Infinity  to  some  minimum  (but 
very  large)  value,  and  then  increases  again.  Thus  tho  neutral  stability 
curve  in  the  a  -  Re  plane  for  aRe  >  >  1  is  a  closed,  isolated  loop.'  Detailed 
studies  bear  out  this  unexpected  conclusion  (Section  -I). 


Since  the  absorption  (or  production)  of  energy  near  the  critical 


layer  becomes  progressively  less  important  with  increasing  M  ,  viscous 
dissipation  must  become  more  important,  in  order  to  counterbalance  the 
production  of  disturbance  energy  near  the  surface.  But  viscous  dissipation 
**(l/aRo),  so  that  the  major  stability  problem  occurs  at  lower  values  of 
«Ro,  whore  the  asymptotic  expansion  procedures  hitherto  employed  may 
become  Inadequate.  For  an  insulated  an  rfacc  K  (Eq.  (13)  approaches  a 
limiting  valuo  of  2  for  >  >  1.  Thus  nRe  at  first  decreases  with 
increasing  Mach  number,  and  then  possibly  "levels  out".  In  that  cass  the 
rapid  decrease  in  Og-v(l/Mc^)  would  eventually  insure  a  rapid  increase 
in  critical  Reynolds  number  with  increasing  Mach  number  at  hypersonic 
speeds.  "V  :.:  '  v 

4.  Eigenvalue  Problem  for  Neutral  Disturbances 

Because  of  the  increasing  importance  of  viscous  dissipation  at  high 
Mach  numbers  and  the  increasing  relative  magnitude  of  the  temperature 
fluctuations,  one  is  justifiably  suspicious  of  ordering  procedures  based  on 
the  well-known  parameter  £  l/(aRe)^.  In  fact  the  full  linearized  disturbanc 

equations  may  be  required.  Nevertheless,  in  order  to  probe  further  into  the 
nature  of  the  high  Mach  number  stability  problem,  we  shall  provisionally 
retain  Prandtl's  "splitting"  of  the  solutions  into  "viscous"  and  "inviscid". 
However,  we  keep  terms  in  the  disturbance  equations  associated  with 
normal  gradients  in  fluctuating  viscosity  and  conductivity,  and  with  viscous 
dissipation,  hitherto  neglected.  For  example,  in  the  energy  equation. the 


he. it  conduction  term  (0/0y+)(  k*  )  glvos  rlso  to  throo  terms, 

Involving  successively  tho  second  and  fir  at  normal  derivative*  of  Q,  and 
0  it  self.  : 

Order  Term 

1  ^i/r)  0"  (the  so-called  "loading  tarm" 

in  the  older  analyses)  ’ 

£  J2/r)|i^*/dT)  O'  T*  (one  contribution  cornea  from  tha 

gradient  In  fluctuating  conductivity, 

~  the  other  from  the  gradient  in 

temperature  fluctuation) 

I  (0/r)(dji/dT)  T’1  (fluctuating  conductivity) 

In  addition  there  in  a  term  -  arising  front  longitudinal  heat  conduction 

(0/Ox*}(k*  );  thl*  term  is  of  order  n*  compared  to  O'*.  If  term* 

of  order  arc  to  be  included,  then  to  be  consistent  one  should  also  Include 
the  a*  £2  terms  arising  from  streamwiso  gradients.  In  the  present  analysis 
only  terms  of  order  1  and  £  arc  retained;  thus,  thld  study  amounts  to  a 
"first-order"  investigation  of  the  high  Mach  number  effects.  This  procedure 
also  has  tha  important  advantage  that  all  viscous  solutions  are  functions  only 
of  the  single  parameter  «Re  (for  a  given  flow),  whllo  thc.lnviscld  solutions 
nre  functions  only  of  a^. 

■  .To  this  approximation  tho  energy  equation  for  the  viscous  disturbances 
is  as  follows: 

Q"  +  (2y^t)(tVu/dT)  TV  *  2(r-l)crMc2  w1!*  -  l-?RPj(w:c)  ft  a  ZsS^lL  0  (23) 

In  the  Dunn- Lin®  analysis  Only  tho  first  and  last  terms  on  the  left-hand  side 
of  Eq.  (23)  appear,  and  the  term  containing  the  normal  velocity  fluctuation 
Is  also  dropped  because  0  is  supposed  to  be  of  order  £  .  Similarly  In  the 
present  study  the  equation  of  continuity  for  the  viscous  disturbances  is 

0'  +  i  f  -  (T'/T)  0  .  My  ,c).  0  a  0  ,  (24) 

and  the  only  term  dropped  is  the  pressure  fluctuation  term,  which  la  of 
order  g  In  the  Dunn- Lin^  analysis  the  0  term  is  also  dropped.  In 
the  Lees- Lin  *  study  this  term  and  the  temperature  fluctuation  term  do 
not  appear,  because  the  analysis  is  limited  essentially  to  cases  for  which 
w  -  c  =  0(  £  ).  In  other  words  the  viscous  disturbances  are  essentially 
incompressible  in  the  Lees- Lin  analysis,  and  the  viscous  temperature 
fluctuations  arc  irrelevant,  so  far  as  tho  eigenvalue  problem  is  concerned. 


In  tho  present  study  the  two  momentum  equations  for  the  viscous 
disturbances  including  terms  of  order  I  and  £  reduce  to  the  following 
single  equation: 

#nii2  d  #,.  ,  i  d  mi  i  l  alto  (w-cj  f  T*  f  _  >i  .  ^  #1  _  n  i 


f»"4  |  T*  f"  -  V -£®J  >  0  .(25) 

In  References  3-5  only  the  terms  containing  fM*  and  f*  are  retained. 

Evidently  tho  three  equations  [Eqs.  (23)  -  (25)]  are  closely  coupled 
and  must  bo  Integrated  simultaneously.  Two  linearly  Independent  seta  of 
solutions  are  distinguished  by  their  bchavlor^n  the  external  flow,  where 
w  and  T  — *■  1.  One  sot  Is  of  tho  form  f  <v  e“  , *  ^‘cl  ^  ,  0  a  0  as 

y  — — ►  oo  »  the  other  set  is  of  the  form  0  **  e"  c)  ^  ,  f  ■  0.  * 

[Tho  third  set  corresponding  to  f*  —^0  as  y—*«»co  corresponds  to  the 
lnviscld  solutions.]  In  order  to  satisfy  the  outer  boundary  conditions  only 
the  negative  exponents  arc  retained.  The  required  numerical  Integrations 
for  M  a  2.2,  3.2,  and  5.6  ware  performed  on  tlia  Data  Iron  205  of  the  Caltech 
Computing  Center  using  a  Rungc-Kutta  Integration  method.  ** 

The  boundary  conditions  at  the  plate  surface  are  that  f,  0,  and  0  : 
all  vanish.  These  conditions  lead  to  tho  complex  secular  equation  for  the 
eigenvalues,  which  after  some  manipulation,  can  bo  written  in  tho  form 

CW  -  l^w'H1*^  (26) 

where  Gw  ■  (ir'/a^ir)w  and  ^  is  an  algebraic  function  of  tho  values  of  the 
viscous  solutions  at  y  a  0.  *  Thus  Gw  ■  Gw  (a;  c,  Me^)  and  ^  ^  (oRej  c,  M#^). 

For  every  value  of  c  the  real  and  imaginary  parts  of  G(ww'/c)  and  (I  -  ^  ) 
are  plotted  in  the  complex  plane,  and  the  intersections  determine  the 
eigenvalues.  . 

Suppose  we  consider  first  the  simpler  situation  at  »  2.  2  and 
limit  the  discussion  to  neutral  disturbances.  Now  -  (ww'/c)(Cq)w  deponds 

mainly  on  vQ(c)  |wc/*w  ,  and  is  practically  independent  of  a.  In  fact. 

It  is  closely  related  to  the  jump  in  Reynolds  stress  across  the  critical 
layer  (Section  3).  As  shown  in  the  sketch  on  page  19,  there  are  two  solutions 

for  each  value  of  -  {ww'/c)(G./w)  >  0.  At  c  =  cfl  this  quantity  vanishes, 


and  one  obtains  the  inviscid  disturbances  a  =  a  ,  aRo 


<x>  ;  there  is  also 


*  For  details  the  reader  is  referred  to  Reference  6. 

**  The  authors  are  grateful  to  Mr.  Kenneth  Lock  for  programming 
and  performing  all  high-speed  digital  computer  calculations. 


lower  branch"  of  the  stability  diagram  for  c  =  c#.  A*  Cy 
e  c a  ,  the  quantity  -  (ww,/c)(Gi)w  Increases  monotonically, 

toward  each  other  on  the  upper  and 

is  reached  for  which  only  one 

MAX 

When  c  <  c 


and  the  eigenvalues  obtained  move 
lower  branches.  Finally  a  value  c  =  c 

solutionis  obtained,  and  beyond  which  no  solutions  exist 
only  one  solution  is  obtained  and  the  rest  of  the  lower  branch  is  traced  out. 
Along  this  branch  Qg  — 0  and  c  — >-  1  -  ( l/Me)  as  Re  — «>  ,  but  aRe  is 
finite. 

In  Figure  5  the  results  of  the  numerical  calculations  are  compared 
-with  the  Laufer ^-Vrebalovich  data  and  the  Lees-Lin3'  4  theory.  Along  the 
upper  branch  (aRe  large)  the  present  calculation  agrees  more  closely  with 


tho  data  {and  will*  the  Dmm-l.ln  calculation)  than  tho  Lecs-Lln  calculation. 
Along  the  lower  branch,  however^  the  values  of  aRo  obtained  by  all  throe 
methods  are  somewhat  too  low,  but  tho  present  method  gives  values  of 
alto  about  30  per  cent  higher  than  tho  Dunn-Lln  method.  Thus  tho 
Inclusion  of  the  terms  of  order  £  seems  to  be  a  step  In  the  right  direction. 
At  one  calculated  point  on  the  lower  branch,  whore  a  0.  030  and 
Re^  ■  89,  one  finds  that  £  a  0.  5d,  so  that  even  more  terms  aro  probably 
required. 

Figure  6  shows  the  calculated  amplitude  distributions  for  one 
eigen-oscillation  on  the  upper  branch,  and  Figure  7  shows  the  mass  flow 
and  total  temperature  fluctuations.  Tho  agreement  Is  surprisingly  good, 
and  confirms  our  belief  that  tho  laminar  stability  problem  at  Mach  numbers 
below  about  2.5  Is  now  fairly  well  understood. 
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At  M0  ■  5.  6  the  situation  is  entirely  different,  at  indicated  already 
In  the  discussion  of  tha  energy  balance  [Eq.  (22),  Section  3]  •  Again 
-  (w  '/cj{G,)  vanishes  at  c  ■  c#  ,  yielding  two  solutions,  one  corresponding 
to  a  ■  ag  and  alts  ■■■>«  «  ,  and  the  other  to  a  finito  value  of  aRo.  As 
Increases  beyond  c  the  quantity  -  (w^'/c)  (G«)  Increases,  but  remains 
very  small,  because  of  the  factor  j  ec/*w|  .  The  upper  Isolated  loop  in 

V** -  Reg  plane  represents  the  solutions  for  uR«  >  >  1  in  tho  region 

cg  *  c  -  1.  The  corresponding  second  solutions  fur  er  *  Cj  ,  c^  ,  Cj  trace 
out  a  new  “upper  loop",  along  which  Op  apparently  0  as  Re  eo  and 
aRc  — finite  value.  Presumably  the  single  solutions  obtained  for 
1  -  (1/Mc)  <  c  <  c  --trace  out  the  rest  of  the  lower  loop. 

The  two  loop  behavior  is  shown  in  Figure  8  and  compared  With 
Dcmetriadcs'  data2.  Since  the  minimum  value  of  Rcfl  on  the  upper  loop 
is  about  10  this  loop  is  of  little  practical  significance.  The  portion  of 
the  tower  loop  drawn  in  an  unbroken  curve  has  about  the  same  shape  as 
Dcmetriadcs'  data,  but  tho  points  lie  about  an  order  of  magnitude  lower 
in  R«p.  However,  tho  calculated  points  are  themselves  about  one  order 
of  magnitude  higher  than  those  obtained  using  cither  tho  Lees-Lln  or 
Dunn-Lln  methods.  Tho  value  of  £  at  the  test  point  indicated  in  Figure  8 
ts  about  2.0,  which  is  some  indication  of  the  inadequacy  of  the  ordering 
procedure  based  on  £  . ;  '■ 

Along  tho  dotted  portion  of  the  lower  loop  in  Figure  8,  c  — l, 
and  ~ — (1  -  c)  may  be  not  be  sufficiently  large  to  justify  the  procedure  of 

splitting  the  solutions  into  "invlscld"  and  "viscous" ^See  Section  2.). 

Numerical  results  were  also  obtained  for  the  "transitional"  case 
Mq  ■  3.  2,  but  these  results  arc  not  yet  well  enough  understood  to  be 
discussed  hero. 

5.  Conclusions  and  Future  Work 

The  present  study  of  the  stability  of  th}  compressible  laminar 
boundary  layer  shows  that  the  relative  Importance  of  the  various  physical  -  - 
mechanisms  governing  stability  changes  drastically  at  high  supersonic 
Mach  numbers. 

1.  Instead  of  being  nearly  constant  across  the  boundary  layer,  the 
amplitude  of  the  inviscid  pressure  fluctuations  decreases  markedly  with 


distance  from  the  plate  surface  at  Mach  numbers  greater  than  3.  Bocauso 
of  this  behavior  the  rate  of  absorption  (or  production)  of  disturbance  energy 
near  the  critical  layer  is  greatly  reduced,  as  compared  with  subsonic  or 
slightly  supersonic  flows. 

2.  At  the  same  time  the  rata  of  production  of  disturbance  energy 
near  the  surface  caused  by  the  viscous  phase  shifts  Increases  with  Mach 
number. 

3.  Viscous  dissipation  becomes  extremely  Important  at  high  Mach 
number,  since  it  must  compensate  for  the  effects  mentioned  In  (1)  and 
(2).  This  phenomenon  is  nlso  foreshadowed  by  the  increase  in  the  relative 
magnitude  of  the  temperature  fluctuations.  Accordingly,  terms  in  the 
equations  of  motion  Involving  gradients  of  viscosity  or  conductivity 
fluctuation,  or  viscous  dissipation,  which  are  neglected  In  the  older  analyses, 
must  bo  included  at  high  Mach  numbers. 

4.  For  free  stream  Mach  numbers  of  2.  2,  and  below,  only  a  single 
stability  loop  In  the  a  -  Re  diagram  is  obtained.  Calculated  neutral 
stability  characteristics  and  disturbance  amplitude  distributions  at 

■  2.2  are  in  good  agreement  with  the  Laufer*-Vrobalovlch  data. 

5.  At  ■  5.  6  two  distinct  stability  loops  are  obtained,  but 

the  minimum  Reynolds  number  (Rc^)  for  the  upper  loop  is  so  high  ( X  10®) 
that  it  does  not  havo  much  practical  significance.  The  other  loop  Is 
qualitatively  similar  to  the  experimental  results  of  Dcmetrladcs^  at 

a  5.6.  However,  the  calculated  Reynolds  numbers  are  still  an  ordor  of 
magnitude  lower  than  the  experimental  values,  although  they  arc  in  turn 

an  order  of  magnitude  larger  than  tho  values  obtained  from  the  Leea-Lln  * 

....  5  •••  .  . 

or  Dunn* Lin  methods. 

6.  At  Mach  numbers  around  3.  5  one  obtains  a  transitional  stability 

diagram  between  the  "almost  Incompressible"  behavior  for  M  =  2.  5  and 

>  e 

the  hypersonic  behavior  for  Mg  =  5.0.  This  regime  requires  additional 
theoretical  study. 

7.  The  structure  and  solutions  of  the  linearized  disturbance  equations 
must  be  carefully  examined  for  the  case  c  — *-  1.  In  addition  there  la  some 
question  concerning  the  existence  of  multiple  eigenvalues  of  the  wave  number 
for  a  neutral  inviscid  disturbance  when  the  relative  velocity  between  the 
wave  and  the  plate  surface  is  supersonic. 
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8.  Asymptotic  methods  utilized  In  all  boundary  layer  stability 
analyses,  based  on  a  ••small"  parameter  of  the  form  £  *  (aRe)"^  ,  aro  no 
longer  ndequalo  at  high  Much  numbers.  In  fact  the  procedure  of  splitting 

the  solutions  Into  "viscous"  (rapidly-varying),  add  invlscld  ("slowly-varying"), 
Is  no  longer  justified.  It  is  suggested  that  the  complete,  linoartsed 
disturbance  equations  should  be  integrated  by  methods  similar  to  those 
developed  in  the  present  study  for  the  separate  viscous  and  Invlscld  solutions. 

9.  Some  evidence  exists  that  aRe  at  first  decreases  with  increasing 
Mach  number,  and  then  approaches  a  constant  as  viscous  dissipation 
builds  up.  Since  the  wave  number  behaves  asymptotically  like  (M#*)*^  , 
the  minimum  critical  Reynolds  number  Is  likely  to  Increase  sharply  for 
hypersonic  speeds. 
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FIG.  I  -  PROPAGATION  VELOCITY  OF  NEUTRAL  INVISCID 
DISTURBANCES  FOR  FLAT  PLATE  BOUNDARY 
LAYERS 
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Ratio  Of  Critical  Point  To  Wall  Pressure  Fluctuation  Amplitudes,  Itc/tw) 


FIG.  3 -PRESSURE  FLUCTUATION  AMPLITUDE  AT  CRITICAL  POINT, 
NEUTRAL  INVISCIO  OSCILLATION 
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FIG.  4  -AMPLITUDE  DISTRIBUTIONS  FOR  NEUTRAL  OSCILLATION 

M#  »  2 .2 


FIG.  7-  MASS  FLOW  AND  TOTAL  TEMPERATURE  DISTRIBUTIONS  FOR  NEUTRAL 
OSCILLATION,  Ma  =  2.2 
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FIG.  8  -  NEUTRAL  STABILITY  CHARACTERISTICS 


